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ABSTRACT 

A  similarity  is  found  betwee^^  the  behavior  of  certain  waves 
in  a  plasma  in  a  magnetic  field  and  gravity  waves  on  the  sxirface 
of  water  of  finite  depthj  in  both  cases,  the  particular  motion 
studied  in  this  report  is  caused  by  an  uniformly  moving  piston, 
started  impulsively  from  resto   This  similarity  in  behavior 
develops  asymptotically  in  time  away  from  the  initial  state,  both 
in  the  linearized  description  of  the  motion  a'nd  in  a  more  general 
non-linear  small-pertxirbation  description  which  includes  the 
linear  asymptotic  motion.   The  approximate  non-linear  theory 
yields  the  third  order  differential  equation  2w^  +  3'ww^  +  'viret:  -   0 
which  can  be  interpreted  as  describing  a  reversible  dispersion 
process,  in  contrast  to  Burgers'  well-known  second  order 
differential  equation  which  describes  an  irreversible  diffusion 
process o 


-  2  - 


NYO-9082 
TABLE  OP  CONTENTS 

Page 

Section 

Xo  XriX/X^OQ,  U.L>  UXOrieoooA    coevoo    ooo*    •o**oooo«    oooooooo    oooeoeeoeoo  /l 

^  O  -^  O  "in  LlX  S.UXOrioooooooooooooo»oooo«o»oeoeoooooooo**oaooeeoo  ~^ 

Jt  O  ijin©9.r      iipprOXXliiS  b  ion  o   <»e«oe«o»««o*e**e«QOOoooa«eeo«oeeoo  J_(J 

[|.o      Non-Linear  Approximat  iono  oo..oeoo«o<»..oooo oooeooo        12 

flppenCl  XX      Xooooooooooooeoooooaoo«a«*o«»««oo«ooeeoo»ooeooo*oeoo  Xy 

ri,ppenCL  XX       XXoo«eo*oooo«ooeooeoo«»oooO9OOoo«o0oo«ttOOoeoo*oeoooO  ^O 

rLCXeX^wnC9Soooo«oooooo9«oeoooeoeo«o«o»eoo»ooooooAooeos«oa«oooo  i^y 


3  - 


NYO-9082 

SIMILARITY  IN  THE  ASYMPTOTIC  BEHAVIOR  OP  COLLISION-FREE 
HYDROMAGNETIC  WAVES  AND  WATER  WAVES 

1.   INTRODUCTION 

In  the  main  body  of  this  report,  we  study  the  motion  of  a  collision- 
free  plasma  in  a  magnetic  field  caused  by  a  piston,  started  impulsively 
from  rest,  moving  with  a  small  vmiform  velocityo   In  particular  we  are 
interested  in  the  asymptotic  time-dependent  behavior  near  the  equilibrium 
stateo   Since  the  study  of  the  corresponding  motion  of  water  waves  is  so 
similar,  the  description  is  given  in  the  Appendices,,   An  interesting 
difference,  however,  is  that  the  original  water  wave  motion  is  two^ 
dimensional  and  becomes  quasi-one-dimensional  asymptotically  for  large  time, 

We  consider  an  Idealized  fluid  model"  of  singly-charged  particles 
moving  in  an  electromagnetic  field.   The  main  properties  of  this  model  are? 
1)   Charge  neutrality  --  the  basic  property  of  a  plasma, 
11)   No  collisions  --  the  charged  particles  interact  only 
through  the  self-consistent  electromagnetic  field,  and 
ill)  Zero-temperature  plasma  --  completely-ordered  motion,  in 
which  the  fluid  motion  is  the  same  as  the  particle  motiono 
The  physical  implications  of  these  assumptions  and  a  study,  primarily  of 
equilibrium  solutions,  are  described  by  Adlam  and  Allen""  (1958a),  Davis, 
Lust  and  Schluter  (1958)  and  other So 


Some  properties  and  solutions  of  this  model  have  been  siiramarized  by 
Gardner  et  al  at  the  Geneva  Conference  (1958)o 

""Also  some  related  calculations  describing  behavior  near  the  initial 
state  of  a  pinched  discharge  are  given  in  a  contribution  to  the 
Geneva  Conference  ( 1958b )o 

Eeg«,  Ao  Banos  and  M,  H,  Mlttelman  at  UCRL,  Llvermore,  and 

Mo  Do  Kruskal  at  Princeton  (private  communication) » 
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The  study  of  water  waves  is  an  old  discipline  aid  the  differential 
equations  governing  their  motion  have  been  known  for  many  years. 
Actually  however,  the  difficult  mathematical  analysis  of  basic  non- 
linear free-boundary  type  problem  is  still  far  from  a  fully-developed 
state,  e.go  see  the  treatise  by  Stoker  (1957)o   The  formulation  for 
water  of  finite  depth  is  summarized  in  Appendix  II  where  the  non-linear 
approximation  is  carried  outo 

The  motion  which  we  study  of  the  plasma  aid  water  waves  may  be 
only  roughly  similar  near  the  initial  state  J  and  the  analogy  between 
the  two  phenomena  develops  only  after  some  time  has  passed  and  the 
motion  is  aoproachlng  the  equilibrium  state.   This  similarity  in 
behavior  exists  asymptotically  in  the  linearized  description  (  ^  3  and 
Appendix  I),  as  well  as  in  a  time-dependent  non-linear  approximation 
(  S   4.  and  Appendix  II)  which  includes  both  the  linearized  asymptotic 
motion  and  the  final  equilibrium  state.   The  range  of  validity  of 
linearized  theory  and  the  approximate  non-linear  theory  is  indicated 
graphically: 

i^ .  _„ ___  Full  non-linear  problem 

<, — ■ Linear  theory  > 


-Approximate  non-linear  theory- 


Initial 
State 


Equilibrium 
State 


->  t 


2o   FORMULATION 

As  a  consequence  of  the  basic  property  of  charge  neutrality,  the 
motion  of  a  cold  plasma  consisting  of  singly-charged  particles  can  be 
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described  by  a  one-fluid  theory^  Yne    state  of  the  plasma  is 

characterized  by  a  system  of  five  first-order  partial  differential 

equations  for  the  variables  (n,  u,  v,  B,  E)  where 

n  =  n_^  =  a_ ,  number  density  of  ions  (or  electrons) 

u  =  u_|_  =  u_,  x-coraponent  of  the  ion  (or  electron)  velocity 

V  =  (v_^-v_)j,  difference  between  the  y-component  of  ion  and 

electron  velocities 

Bj  z-coraponent  of  the  magnetic  field 

E,  y-component  of  the  electric  field 

(x  j 
and,  with  the  exception  of  E    ^  x-component  of  the  electric  field 

(charge  separation  field),  all  other  comnonents  are  zeroo   These  six 
quantities  depend  on  the  single  space  coordinate  x  and  the  time  tc   For 
convenience  we  consider  these  quantities  to  be  normalized  in  the  follow- 
ing way; 

ns   n 
o 

ut   A  =  [B^/^n  (m^  +  "i J  ]  "^^^ 


o     o'  ^  o   +    - 

v°   a  =  [B  (m^  +  ra  )/un  m.m  j  ' 
O  0+     -o  +  - 


B:   B 
o 

E:   E  =  A  B 
c    o  o 


E^^^  e(^^=  a  B 
o     o  o 


x;  X  =  [m.m  /ue^  n  (m, +m  ) ]  ' 
o      +  —  "^+0+   — 


where  \i   is  the  magnetic  permeability  of  space,  e_^  and  m^  are  the  charge 
and  mass  of  the  ion  (minus  sign  for  electrons),  and  the  subscript 
(zero)  refers  to  a  constant  initial  or  boundary  state.   Then  the  re- 
sulting normalized  equations  are  one  conservation  eauatlon  of  number 
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density  and   two   of  moTnentum  and  Maxwell's   equations    (two): 

n^   +    (nu)      =  0  (1) 

t  X 

u .    +    uu     =   vB  ( 2 ) 

t  X 

v^   +   uv     =  E   -   -uB  (3) 

t  X 

B.    +   E     =   0  ik) 


B      =   -nv  (5) 

X 


where  the  displacement  current  is  neglected  in  (5)  and  lettered  sub- 
scripts indicate  partial  differentiation.  In  this  one-fluid  theory, 
the  x-component  of  the  electric  field  is 

/  \    III,  ~  in 
e(^)  =  (_!■  ^)vB  (6) 

rri,+  m 
+    — 

Also  by  eliminating  the  y-component  of  the  electric  field  E  between 
(3)  and  {I4.)  we  get  a  simplification.   Then 


(v  +  B).  +  [u(v^  +  B)]^  =  0  (7) 


Comparing  (7)  with  (l),  (v  +  B)  satisfies  the  same  equation  as  n, 
Thus  we  can  set 


V  =  n  -  B  (8) 

X 
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if  this  is  true  initiallyo   Now  we  have  a  reduced,  but  equivalent^ 

system  of  four  equations  (l),  (2),  (5)  and  (8)  for  (n,u,v,B)  in  place 

of  the  original  five  and  (3)  may  be  regarded  as  the  defining  equation 

for  Eo   We  note  that  the  particle  path  and  the  plane  t  =  constant  are 

the  only  set  of  real  characteristics  (actually,  a  double  set)  for  this 

differential  equation  system^  and  neither  A  nor  a  are  characteristic 

''     *  o      o 

velocities  in  this  sense,, 

For  later  use  {^  l\.)    we  look  for  equilibrium  solutions  of  the 

reduced  system  (1),  (2),  (5)  aid  (8)5  noting  that  the  equations  are 

Gallilean  invariant,  ioeo,  invariant  under  the  transformations 

fx  ■— *■  (x=u  t)  o  t  — »>t,  u  — *•  (u-u  )[  where  u  is  constanto   Setting 
(_  o  o_^        o  ° 

a 

masSj,  momentum  and  energy  follow  by  inspection 


-^  =  0  in  (1)  and  (2),  combined  with  (5)  and  (8),  three  Integrals  of 


nu  =  u  (9) 

o 

2 

nu  +  ~-=   xx/   +   I  (10) 

O     <d     o     ^ 

^  2 
~{u^   +  v^)  +  B  =  -|-  +  1  (11) 

Combining  (10)  and  (11)  yields  the  "phase  relation" 

v2=  (B.l)2[l  -  ^^ii^  ]  (12) 

Eliminating  n,  u,  and  v  by  (5),  (9),  (10)  and  (12),  B(x)  can  be 
integrated  explicitly  (cfo  references  in  ^  1 ) o   We  obtain  the  solution 
in  the  limiting  case  of  a  small  disturbance  by  the  transformation 


8  = 


?  =  e^/^x  and  letting  n  =  1  +  en^^^^),  u  =  1  +  eu^-^^(?),  v  =  e^/^v^  "^^  (4) , 
B  =  1  +  eB    (<?)  and  u  =  1  +  ea,  where  a  =  const,  to  be  determined  and 
e  =  (B    -1)«  1  is  the  magnetic  field  amplitude.   Then  (5),  (9),  (10), 
and  (12)  become 


f^=-v'^'  (13) 


U<1'  +  B<1'  =  a  (15) 


(v<l')2=  (B'^')2(2a  -B<1))  (16) 


Eliminating  v^"""^  between  (13)  and  (16),  B^-"-^  satisfies 


||!H=  -B(l)(2a-  B^1))V2  (,^j 


which  yields  the  equilibrium  solution 

B^^^  =  2a  sech^(a?)  (l8) 


We  choose  a  =  l/2  so  that  B    =1,  consistent  with  our  definition  of 

max     ' 

ta       (l8)  is  a  pulse-type  solution  similar  to  the  small-amplitude 
solitary  water  wave,  e.go  see  Keller  (I9I4-8),   The  bouidary  conditions 
chosen  in  (9),  (10)  and  (11)  have  eliminated  some  additional  equilibrltoui 
solutions  (periodic)  which  can  be  recovered  from  the  time-dependent 
non-linear  approximation  derived  in  S   Ij-o 
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3,   LINEAR  APPROXIMATION"" 

To  describe  the  motion  caused  by  small  initial  distiirbances  we 
linearize  the  reduced  plasma  equations  (1),  (2),  (5)  and  (8)  by 
perturbing  on  the  rest  state  (n,u,v,B)  =  (l,0,0,l)o   The  resulting 
equations  in  the  first  approximation  are 


njl)  +  u  ^^^  =  0  (19) 


u^^l)  =  v^l)  (20) 

B  (1)  =  -v^l)  ,  (21) 

X 


V  (1)  =  n^l^.  b(1)  (22) 

X 


By  elimination  we  obtain  a  single  fourth  order  differential  equation 

(1) 
for  the  x-component  of  the  velocity  u    ; 


Also  n^"''  ,  B  ^^  and  v^"*"^  satisfy  the  same  equation  (23).   We  now 
consider  the  plasma  motion  for  x  >  0  caused  by  an  electrically- 
neutral  piston'"""'  (nominally  at  x  =  0)  started  impulsively  from  rest 


"      The  results  in  this  section  were  presented  by  Gardner  (1958)  at  the 
Controlled  Thermonuclear  Conference  in  V/ash«  D,  C.  and  have  been 
summarized  in  the  Proceedings  of  the  Geneva  Conference  by  Card  ler  et  al 
(1958)c 

This  idealized  problem  is  considered  for  convenience.   We  could  just 
as  easily  take  the  same  differeitial  equation  (23)  satisfied  by  B^-"-^ 
and  consider  a  magnetically-driven  pistono 
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and  moviig  with  uniform  velocity  U,  in  the  positive  x  direction.   Then 
the  initial  conditiois  are 

u^^^x,0)  =  u^^^^x,0)  =  0  (21+) 

and  the  boundary  conditions  for  t  >  0  is 

^ ^^-^  =   H(t)  (25) 

1 

where  H(t)  is  the  unit  step  function  (=  0  for  x  <  0  and  =  1  for 
X  >  0)o   Equation  (23)  with  the  conditions  (2L|.)  and  (25)  can  be  solved 
by  applying  the  Laplace  transform  with  respect  to  the  time  t.   The 
solution  is 


"  U^"'"  -5iT   I   f  •'^P  [ts(l-a(l.s2)-l/2l]      (26) 


.(1) 

T       2Tii 

T 

where  a  =  x  /t  and  the  contour  /^  goes  from-ioDto  +  i  qd  in  the  right 
half  of  the  complex  s-planeo   To  carry  out  the  inverse  Laplace 
transform  (26)  apoears  difficult^  but  the  asymptotic  behavior  for  large 
t  can  be  determined  by  the  method  of  steepest  descent.   Then  the 
asymptotic  form  of  the  velocity  distribution  for  large  t  (and  x)  is 


( 1 )  /   ,  ^     roo 
P 


where  p  =  (2/3)"'-^^(x-t  )/t-^/^  and 
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Ai(p)  =  ^ 


poo        ,3 

dX  cos(-^  +  SX)  (28) 


is  the  Airy  functiono   The  same  asymptotic  behavior  is  fouad  for  the 
linearized  motion  of  water  waves  in  Appendix  I,  equation  l-17o   Equation 
(27)  describes  a  shock- type  velocity  distribution  with  damped  oscilla- 
tions behind  the  "shock"  front  (p  <  0)o   The  velocity  decreases 
exponeatially  ahead  of  the  front  P  >  0,  and  for  negative  (3  the  velocity 
oscillates  about  and  approaches  the  piston  velocity  U,  as  p  —^    -00  0 

_0  /l 

Also  the  amplitude  of  the  oscillation  decreases  like  p"-^^' ^  as  p  -*•  -oo  c 

However  this  wave  is  not  a  shock  in  the  usual  sense  since  it  is  non- 

1/3 
stationary,  widening  with  time  like  t  '  °  thus  (2?)  describes  a  slowly 

dispersing  wave© 


i^o   TON-LINEAR  APPROXIMATION 

So  far  the  only  apparent  relationship  between  the  linearized 
asymptotic  behavior  of  plasma  motion  (equation  (27)  in  ^   3)  and  the 

pulse-like  equilibrium  solution  (equation  (I8)  in  o  2)  is  their 

2  1/2 

propagation  speed  A   =  [B  /un  (m,  +  m  ) ]  ^   as  e  —^    0©   In  this 
^"  o     00+    - 

section  we  seek  that  approximate  time-dependent  description  which  is 
valid  for  large  time,  ioeo,  valid  under  the  conditions  that  it 
1)  includes  the  equilibrium  solution  (I8)  and  2)  is  in  the  same 

asymptotic  direction  in  the  x-t  plane  as  the  linearized  asymptotic 

1/3 
solution  (27),  that  is,  the  direction  defined  by  p  =  (x-t)/t  ■ 

consto   It  is  clear  from  (17)  and  ( I8 )  that  a  non-linear  approxima- 
tion is  requlredo   A  comparable  formulation  has  been  found  in  the 
approximate  description  of  the  motion  of  flood  waves  in  rivers  by 
Morikawa  (1957). 
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To  obtain  the  desired  non-linear  approximation  of  the  plasma 
equations  (1),  (2),  (5)  and  (8),  we  first  apply  the  above  condition 
1)  which  states  that  the  pulse  solution  (l8)  is  to  be  included  in 
our  approximate  formulation.   Thus,  we  make  the  following  scale 
transformations : 

C  =  e'^^^ix-t)  (29) 

and 

V  =  e^/^V.  (30) 

Condition  2)  above  fixes  the  necessary  scale  tra-isformation  of  time: 

l/3 
since  we  want  (x-t)/t  '   to  be  invariant  with  respect  to  the 

coordinate  transformations  (29)  and  (31).   In  addition,  we  make  a 
perturbation  expansion  with  respect  to  e  on  the  rest  state  (n,u,v,B)  = 
(1,0,0,1): 

n  =  1  +  en^^^?,r)  +  e^n^^)^^^^)  +  ^_  (32) 

u  =  eu^^^?,r)  +  e^^^h^sT)   +  ...  (33) 

V=  eV^^^^,r)  +  e^v^2^c,r)   +  ...  (31^.) 

B  =  1  +  tB^^U^,r)  +  E^B^^U^,r)  +  ...  (35) 
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By  the  scale  transformations  (29),  (30)  and  (31),  the  equations 
(1),  (2),  (5)  aad  (8)  become 

er^  +  [n(u-l)]^  =  0  (36) 

e^^   +   (u-l)u^  =0  (37) 

B^  =  -nV  (38) 

eV^  =  n-B  (39) 

Putting  (32  to  35)  in  (36  to  39)  and  equating  the  coefficients  of 
powers  of  e  in  the  usual  way,  we  obtain  the  equations  which  the  lowest 
order  approximation  (n^-*"  ,  u^-'-  ,  V^-"-  ,  B^''-')  satisfies: 


n'l'  -  u'l'  =0  (1^0) 

41*  =  B^l'  =  -v'l'  ((H) 


and 


(n(2'-u<2))^  =  4l'  .  (n'l'u'l))^  (W) 

u<2'  *   V<2'  =  41'  *   u<l'u(l)  -  V<1>b(1>  (W 

b(2)  +  v(2'  =  -nd'v'l'  (U5) 

n<2'  -  B<2'  =  v'l'  (1;6) 

(?)   (2)   (2)   (2) 
The  second  order  approximation  terms  (n^   ,u^   ,V^  ',B^   )  may  be 

eliminated  among  (1^.3  to  1+6)  and  u^-""^  (or  n^-""^  or  B^-^M  satisfies  the 

third  order  non-linear  equation 
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2u<lU   3u<l'u<lU  4U  =   0  (^7) 


(14.7)  describes  the  plasma  motion  traveling  at  the  velocity 

A  =  [B  /un  (m.  +  m    )]   '       In  the  positive  x  direction  near  and  in- 
o     00+    — 

eluding  the  equilibrltjm  state.   An  equation  of  the  same  form  describes 
the  motion  of  water  waves  traveling  at  the  speed  of  small  gravity 
waves  (Appendix  II,  equation  ( 11-32) )• 

In  contrast  to  Burgers'  second  order  equation  which  describes  an 
irreversible  diffusion  process,  (l^.?)  describes  a  dispersion  process 
which  is  reversible  in  the  following  sense:   the  equation  is  invariant 
with  respect  to  change  of  sign  ^  — *■  -B,   and  "^ — >    -to      The  solution 
of  the  initial  value  problem  with  a  step-function  initial  distribution 
u^-'-^(?,0)  =  I  -  H(?)  (H  =  0  for  4  <  0  and  H  =  1  for  ?  >  0)  would  be  of 
some  interest.   But,  the  Hopf  (1950)-  Cole  (1951)  technique,  which 
yields  the  corresponding  solution  of  Burgers'  equation,  appareatly  does 
not  apply  here.   We  plan  to  carry  out  a  numerical  integration* 

'  It  is  easily  shown  how  the  linearized  asymptotic  solution  (2?)  is 
related  to  (1|7)  and  that  the  pulse  solutioi  (I8)  is  included  in  the 
equilibrium  solutions  of  (lj-7).   Neglecting  the  non-linear  term  in 
(1^.7)  we  look  for  similarity  solutions  of 

2uj,^^  +  u^J^  =  0  (1^8) 

of  the  form  u    =7r™f(S)  where  (3  =  bCt"^  and  b,  m  and  n  are  constants* 

Then,  m  is  arbitrary  and  n  =  -l/3  and  f  satisfies  the  ordinary 

differential  equation 

<{. 

Burgers'  equation  is  w.+  w  w  =  v  w   where  v  =  const. 
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t  J  I      J  ,,  , 

f    =  pf  +  3mf  =  0  (l|9) 


1/3 
where  b  =  (2/3)    for  convenienceo   We  consider  three  particular 

values  of  m  =  0,  -l/3,  =2/3: 

Case  1,  m  =  Op   For  this  case  {l\.9)    becomes 


(f')"  -   p(f')  =  0  (50)    ! 


and  f   satisfies  the  Airy  equation  so  that 


f'  =  A^(P)  (51) 


The  linearized  as3rmptotic  solution  (2?)  corresponds  to  this  casoo 
Case  2,  m  =  -l/3«   For  this  case  a  solution  of  (1^9)  is 

f  =  A^(P)  (52) 

and 

^(1)  ^^=1/3  ^^(p)  (53) 

(53)  is  the  asymptotic  solution  obtained  by  Jeffreys  aid  Jeffreys 
(1956)  in  their  study  of  the  motion  of  water  waves. 
Case  3.  m  =  °2/3o   For  this  case  a  solution  of  (L|.9)  is 


f  =  a/(b)  (51;) 

i 
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(53)  is  imbounded  for  p  — >  -co but  this  case  is  of  interest  since 
the  non-linear  equation  (I4.7)  also  has  a  similarity  solution  of  the 
form  u    =T"^/-^f(p)  where  f(p)  satisfies  the  equation 


f   +  ff  -  pf  -2f  =  0  (55) 


It  would  be  of  interest  to  determine  whether  there  is  a  bounded 
solution  of  (55)  over  the  entire  range  of  -00  <  P  <  oo.  Equilibrium 
solutions  of  (1+.?)  are  obtained  by  seeking  solutions  of  the  form 
u^  '  =  P(^-a'C)  where  a  =  const,  to  be  determined.   Since 
(^-ar)  =  z^'^lr.  -   (l+ea)t],  (ea)  is  the  velocity  deviation  of 
the  steady  progressing  wave  from  the  critical  speed  A  •   By  (I4.7), 
F(C  -  af)  satisfies  the  equation 

P'"  +  3Ff'  -  2aF  =  0  (56) 


Multiplying  (56)  by  F  and  integrating  twice  we  obtain 


(p')2  +  p3  _  2aF^  =  c^^P  +  C2  (57) 


where  c,  and  Cp  are  integration  constants.   If  we  impose  the  condition 
that  the  solution  vanishes  at  infinity,  then  p  =  P  =  P  =0  implies 
that  c,  =  Cg  -  0;  and  (57)  becomes 

(p')2  =  p2(2a  -  p)  (58) 
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which,  corresponds  to  the  solitary  wave  equation  ( 1? )  o   If  we  define 
e  as  the  wave  amplitude,  the  a  '^}.'   =   1  implies  that  a  =  l/2o   The 
flow  velocity"  is  supercritical  at  infinity  and  subcritical  at  the 

peako   Periodic  solutions  are  obtained  by  relaxing  the  conditions 

« 
at  infinityo   For  example  P  =  P   =0  implies  that  c-  =  0  and  we  obtain 

for  fixed  amplitude  a  one-parameter  (wave  length)  family  of  periodic 

solutions  (cnoidal  waves)  with  speeds  less  than  or  equal  to  the 

solitary  wave  speed„ 

Pinally,  we  summarize  the  correspondence  between  variables 

describi^.g  the  asymptotic  motion  of  plasma  waves  and  water  waveso 

Plasma  Water 

t   =  X  /A  ^K/s)^"^^ 

o    o'  o  o  "^ 

We  point  out  that  for  plasma  motion  the  y-component  of  velocity 

( 1  ) 
V    is  normalized  with  respect  to  a  different  velocity 

a  =  [B   (m^+  m  )/u,n  m_,m  ]^       than  u    which  is  normalized  with 

respect  to  A  (cfo  6'  2 )  o   Thus  the  actual  y-component  of  velocity 

exceeds  the  x-component  of  velocity  by  a  factor  a  /A     =    (m^m_)    o 

[1  +(m_/m^)]  for  (m^m^^)  »  lo 


*""  Note  that  the  flow  velocity  is  from  left  to  right  with  respect  to 
the  solitary  wave  in  ^  2  and'  in  the  opposite  direction  in  this 
section,  appropriate  for  a  piston  moving  from  left  to  righto 
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APPENDIX  lo   LI^IEAR  APPROXIMATION  OP  MOTION  OP  WATER  WAVES 

To  compare  with  the  plasma  motion  described  in  6   2,  we  consider 
the  motion  caused  by  a  piston  (started  impulsively)  moviig  at  a  small 
uniform  speed  U^  into  water  initially  at  rest  of  uniform  depth  h^o 
The  piston  (nominally  at  x  =0)  moves  in  the  positive  x  direction 
and  the  gravitational  acceleration  g  is  in  the  negative  y  directiono 
Neglecting  viscosityj  the  motion  is  described  by  the  velocity  potential 
j25(xyy,t)„   For  convenience  we  introduce  the  following  normalization: 

x;    h 
o 

t:    (hyg)l/2 

Since  the  boundary  conditions  on  the  water  surface  (cfo  Appendix  II) 
are  linearized,  the  flow  is  nominally  confined  to  move  in  a  semi- 
infinite  strip  x  >  0,  with  the  horizontal  bottom  at  y  =  0  and  the 
free  surface  at  y  =  lo   The  motion  is  governed  by  the  Laolace  equation 
in  the  Interioro 

0      +  j6      =  0  (i-i) 

'^xx   ^yy 

The  initial  conditions  are 

0(x,y,O)  =  ^^(x,y,0)  =  0,  (1-2) 

The  boundary  conditions  are:   at  x  =  0 

^     ~  =  H(t)  (1-3) 


"l 
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where  H(t)  is  the  lonit  step  function^  on  the  bottom  y  =  0 


0^(x,o,t)  =  o; 


(I-It) 


and   at    the  linearized  free    siirface   y  =   1 


S^^tt  ^  ^y  =  0 


(1-5) 


The  solution  of  (I-l),  with  the  conditions  (1-2  to  5),  is  readily 
obtained  by  applying  the  Laplace  transform  with  respect  to  t  and  the 
Foiorier  cosine  transform  with  respect  to  xj  the  resulting  second 
order  ordinary  differential  equation  in  the  remaining  variable  y  is 
easily  solvedo   After  a  simple  inversion  of  the  Laplace  transform. 


the  inverse  Fourier  transform  yields  the  integral  solution 


00 


^{yi,Y,t) 


2U, 


■K 


1 1/2, 


Jo 


T^  cos(^^^)   1-  ^fg^^  cos[^tanh  <)"/"t]y(I-6) 


The  wave  profile  n(x,t)  =  0. (x,l,t)  and  from  (1-6) 


oo 


2IL  ( 


y(x,t)  =  -  ^  '    ^X(K.tanh  )^)^/^cos(i^<:)sin[^;tanh|^)^/^t]       (1-7) 


We  make  an  asymptotic  evaluation  of  (1-7)  for  large  t  (and  x)  by  a 

variant  of  the  method  of  steeoest  desceito   Since  the  integrand  of 

f°°   1 
(1-7)  is  regular  at  )( =  0  and  is  an  even  function  of  /(_,  then  1   =  ^ 

~^-oo 

We   xririte    (1-7)    as    the    s-'om   of   two    sine    Integrals 


J 


OD 


•OO 


where    a  ^'Y^x;,    b   =    (  K.tanh  Y^)         t    an 


cos    a    sin   b   =   -   i[sin(a-b)    -    sin(a+b)] 
1/2. 


(1-8) 


20 


(a+b)    =Xt[a   +    (tanhXy-^/^/X,^/^]  •  (1-9) 


where    a  =  x/to      For  ){«    1    (large   x),    (1-9)    is    approximately 

.2 
{a+h)^tt[a  +    (1-  ^-)]  (I-IO) 

For  lar^e  t  and  a  =  0(1)  (x  >  0),  the  Integral  containing  sin(a-b) 

l/3 
dominates  that  with  sin(a+b)o   By  transforming  ;i;^=  (2/t)  ^  ^X    and 

(a-1)  =  p/(2t^)-^/^,  we  rewrite  (a-b)  in  (I-IO) 
(a-.b);^^(^  +  P) 


and  the  asymptotic  solution  of  (1-7)  for  large  x  and  t  is 

\]      r^  3   _ 

7l(x,t)^2^  i    ^  sln(^  +  p\)  (l-ll) 


J 


CO 


Since  the  derivative  with  respect  to  p  of  the  integral  (I-ll)  is 
the  Airy  function 

QO 

Ai(p)  =^i  d\    cos(^+p^)  (1-12) 


we   rewrite    (I-ll)    in   the    concise    form 

/^oo 
/\    (x,t)    ^  TJ^  Ai(Y)dY  (1-13) 

^  P 

Xv'here 

P   =   2l/^(x-t)/tl/^  (I-II^) 
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The  i/'ave  profile  (I-lji)  Is  of  the  same  form  as  (2?)  in  •.'  2o 
Ccmpai-lng  ( I-II4.)  with  (28),  p  differs  from  (3  by  a  factor  3  -^^ -^ 
(cf  Appendix  II,  equation  (II-32 ))  »   A  more  direct  comparison  with 
(27)  is  found  by  calculating  the  x-component  of  velocity 
u(x,l5t)  =  ^    (x,l,t)  at  the  free  surface  y  =  1,   Prom  (1-6) 

^^g^^>^^  =\j        ^   sin(Xx)  fl-cos[(Xtanh  <)^/2^]j  (I-l5) 


Since 


I    -^  sin(/l^)  =  Tl,  X  >  0, 
J  -co 


(l-l^)  becomes 


^^f/^>^^  =  1.1      %-^  sin(Xx)  cosLXtanh  /J^/^t]  (l.-^^) 

u,         ^  j   X. 

-^  -00 

which  differs  from  the  wave  profile  (1-7)^  but  the  asymptotic  form 
of  (1-16)  is  the  same  as  (1-13)  and 


CO 

u(^,\itAl  ^^  ,  Ai(Y)dY  (1-17) 

^1        J 
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APPENDIX  II o   ASYMPTOTIC  NON-LI ^AR  APPROXIMATION  FOR  MOTION  OP  WATER 
WAVES 

In  two  dimensions  the  system  of  equations  governing  incompressible 
flow  are  the  conservation  equations  of  mass  and  mfementum: 

u^+Vy=0  -  (II-l) 

ff  -  P,  =  0  (II-2) 


ffH.p^.l  =  0  (11-3) 


where 


at       t       X       y 


The  variables  have  been  normalized  as  follows; 


x,y-  h^ 


t^  (hyg)^/2 


pt  pgh^ 


where  h   is  the  depth  for  the  rest  state,  g  is  the  gravitational 
acceleration  and  p  is  the  density  (constant),,   In  addition  to 
(11-1,2,3)  we  impose  the  irrotationality  condition 

V   -  u  =  0  (II-1|) 

X    y 

The  boundary  condition  at  the  horizontal  bottom  is 

v(x,o,t)  =  o;  (II-5) 
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and  at  the  free  surface 

v(x,h,t)  =  I  ^  +  u(x,h,t)  ^  [  h(x,t)  (11=6) 


[^^^-'^ 


and 

p(x,h,t)  =  Oo  (II-7) 

Just  as  for  the  plasma  equations  (l),  (2),  (5)  and  (8),  the  particle 
path  and  the  plane  t  =  constant  are  each  double  characteristics  of 
(II-l,  2,    3)o 

Following  the  same  procedure  by  which  we  derived  the  non-linear 
approximation  of  the  plasma  equations  in  c  l\.,    we  obtain  the  asymptotic 
approximation  for  water  waves  by  1)  including  the  small  amplitude 

solitary  wave  as  an  equilibrium  solution  and  2)  scaliig  the  time  so 

l/3 
that  (x-t.)/t  '   is  invariant p  as  dictated  by  the  linearized  asymptotic 

solution  (I-17)o   Then^the  required  scale  transformations  are  the  same 

as  those  for  the  plasma  equations  ( cf o  (29),  (30)  and  (31)) 

V  =  e-^/^V  (11-10) 

and  the  remaining  variables  yj,u  and  p  are  vmchanged.   The  flow  equations 
(II=l,2p3)  are  transformed  by  (11-8,9,10)  into 

Uj.  +  V  =0  (11=11) 

1^  +  p^  =  0  (11=12) 

e  H  +  Py  +  1  =  0  (11-13) 

-  2k  " 
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where  the  notation  d(  )/dt'  means 

^P-  =  e(  L  +  {n-l)i    ).  +  Vo(  )  (Il-llj.) 


and  the  irrotationality  condition  (ll-l^.)  becomes 

eVg  --  u  =  Oo  (II-15) 

In   addition^    the   boundary  conditions    (ll°5»6j,7)   become 

VlCpO.r)    =  0  (11-16) 

V(^,,h,t)    =      e  ^  +    [u(C,h,r)-=   1]   ^  I    h(^,T)  (11-17) 

and 

vUph.T)    =  0  (11-18) 

There  is  a  close  relationship  between  the  scale  transformations 
(Il=8j,9pl0)  and  the  transformation  of  variables  yielding  the  well-knovm 
long  wave  (or  shallow  water)  approximation  [cfo  Friedrichs  (19I4-8)  or 
Morikawa  (1956)] o   Both  transformations  yield  the  hydrostatic  approxi- 
mation  for  e  -»  Op  ieeo  in  (11-13)  e  ^  «  (p  +  Do   Also  in 
Morikawa" s  version  (differs  from  Friedrichs'  in  the  scaling  of  the 
vertical  velocity  component)  the  irrotionality  condition  (11-15)  states 
that  u  =  0  for  6=0^  ioe,,  the  horizontal  velocity  component  u(4,T) 
is  independent  of  the  vertical  coordinate  yo   The  scale  transformation 

yielding  the  loig  wave  approximation  differs  only  in  the  time  scaling, 

1/2 
being 't^=  e  '  t  instead  of  (II=9)o   That  is,  the  non-linear  approxlma- 

3/2 
tionj,  which  we  describe  in  this  Appendix  ('t'=  e  '  t),  describes  the 
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motion  at  later  time  and  closer  to  the  final  equilibrium  state  than 
the  long  wave  approxlmatlono 

In  addition  to  the  scale  transformation  we  make  a  perturbation 
expansion  (cfo  S   1|)  on  the  rest  state  (hpU,V,p)  =  (1,0,0,  (l-y))o 
Then, 

h  =  1  +  eh^^h^.T)   +  e^^^^Cpt)  +  ooc  (11-19) 

u=  eu^^^CpypT)  +  e^u^2^4py,t)  +  co.  (II-20) 

V=  eV^^^(5,ypr)  +  e^V^2)(^^y^^)  +  ^^^  (11-21) 
P  =  (1-y)  +  e[p^^^?,y,T)  +  h^^h^.^i)]   +  e^lp^^U^.j,^)   + 

+  h^^U^.r)]   +   ooo  (11-22) 

noting  that  at  the  free  surface  y  =  h 

v(?,hpr)  =  ev^i^?,ipr)  +  e^[v^2)(^^3^^^)  +  vi^h^pi,r).h^^^^pr)]+ooo 

(11-23) 


p(C,hpt')  =  ep^^^€,l.T)+  e^[p^^^4plpr)+  p^^^4.1,T)oh^^^CpT)]+  oo, 

(II-2i|) 

Putting  (11-19  to  22)  in  ( 11-11  to  l8)  and  equating  coefficients  of 
like  powers  of  e,  the  lowest  and,  if  desired,  higher  order  approxima- 
tions can  be  calculated  in  a  systematic  wayo   We  summarize  the 
equations  yielding  the  lowest  order  approximation  (h   ,  u   ,  V   ^ 
p^-^^)o   Thenp  u^^U^.T)    is  independent  of  y,  p^-'-^E  0  and 

h^l^  -  u^l^  =  0  (II-25) 

v^l)  -  -y  u^^)  „  (11=26) 
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AlsOp 

u(2)  ^  ^(y2/2)u^l)  +  p(^^^)  (11=27) 

V^2)  =  (y3/6)u^y  ==  yF^  (11-28) 

p(2)  ^  l(l„y2)u^l)  (II-29) 

4^^  =  "[^4M/^^  "  ^^^  ■"  4^^  ""  (h^^^^^^^)^  (11-30) 

(p(2)  ^  t^(2)  .  u^2)^^  ^  ^(^|1)  ^  ^(1)^(1))  (II„31) 

(2)    (2)    (2)    (2) 

The  second  order  terms  (u  \    v    p  p    p  h    )  and  the  integration 

variable  (with  respect  to  y)  F(4sT)  can  be  eliminated  among  (11-27  to 

31)  resulting  in  a  single  equation  for  u    ; 


2u^l^  +3u^^)4^)  ^  l4cC=  °  ^^^  =  ^2^ 


Except  for  the  factor  (l/3)  multiplying  the  third  derivative  term  and 
which  can  be  eliminated  by  renormallzing  the  space  and  time  variableSp 
(11=32)  is  the  same  equation  as  ( I4.7 )  which  describes  the  asymptotic 
plasma  motiono   Aside  from  minor  differences,  (11=32)  was  derived  many 
years  ago  by  Korteweg  and  de  Vries"(l895)  following  an  approximation 
procedure  of  Rayleigh"  our  procedure  is  more  direct  and  indicates  more 
clearly  the  asymptotic  nature  of  ( 11=32  )o   Like  ( L|-7 )  for  plasma  wavesp 

the  third  order  non-linear  partial  differential  equation  (11=32) 

1/2 

describes  the  motion  of  water  waves  traveling  with  velocity  (gh  )  ' 


Keller's  (19I4.8)  results  are  also  a  rederivation  of  the  small-amplitude 
equilibrium  solutions  by  scale  transformation  with  respect  to  a  small 
parameter e 

=  27  = 


in  the  positive  x  direction  near  and  including  the  equilibrium  state, 
A  related  fourth  order  differential  equation  which  approximately 
describes  waves  traveling  in  both  positive  and  negative  x  directions 
due  to  Boussinesq  has  been  rederived  recently  by  Ursell  (1953) o 
(11-32)  can  be  derived  directly  from  this  fourth  order  equation  by 
applying  the  scale  transformations  (11-8,9)  and  letting  e  — >  Oo 
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